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ABSTRACT. In this work, we present reproducing kernel method to solve non- 
linear system of differential equations. Several numerical results are given to 
show the efficiency of the method. 
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1. INTRODUCTION 


System of differential equations and dynamical system have been widely used in 
applications such as mathematical model, statistics, biological application, popula- 
tions etc. Consider a system of differential equations 


w(t) = f(t, a(t), y(t)); 


y'(t) = g(t, x(t), y(t)); 
subject to 2(0) = 0 = y(0), where f and g are continues in t, x, and y. 


(1.1) 


2. PRELIMINARY 


This section is devoted to some basic definitions, lemmas, and theorems which 
will be used in later. 


Definition 2.1 ( [1]). Let H be a Reproducing kernel Hilbert space of continuous 
real-valued functions f : S —> R. Then there exists k: S x S > R in H such that 
(f(-),kC, 8)) = f(s), for all f €H and all se S. 


Definition 2.2 ( [1]). Let W3[0,7] be the space of all functions u : [0,7] > R 
such that u is absolutely continuous on [0,7] and its derivative belongs to L?[0, T]. 


Definition 2.3. Let W2[0,T] be the space of all functions v : [0,7] > R such that 
v' is absolutely continuous on [0, T] and v” belongs to L?[0,7']. 


Definition 2.4 ( [2]). The space W4[0, T] equipped with the inner product 
(u,e}wgqocry =O) + f (odr, 
[0,T] 
is a RKHS with reproducing kernel function q(-,-) given by 
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(2.1) q(T) = ane eee) 


Email addresses: jshm97@mst.edu 
2010 Mathematics Subject Classification. 46E22,47B32, 47B38. 
Key words and phrases. reproducing kernel Hilbert space, dynamical system. 


1 


2 JABAR S. HASSAN 
Definition 2.5. [1,2] The space W2[0, T] equipped with the inner product 
(ue) gtocry = w(O)0(0) +w (OO) + fw" (ryu"(r)ar, 
[0,7] 


is a RKHS with reproducing kernel function r(-,-) given by 


(2.2) later 


Definition 2.6. [2] Let W’[0,T] = W2[0,T] 
vector functions fi : [0,7] x [0,7] > R where 
for j = 1,2. 


® W3[0,T] be the space of all real 
h = [hi,h2]" and hy, ho € WE[0,T] 


3. REPRODUCING KERNEL METHOD 


In this section, the general structure of representation of the exact solutions will 
be discussed to problem (4.1). 

We define operators F,G : W3[0,T] + W3[0,7], for each x,y € W3[0,T] such 
that 


| Fe()=2"() 

(3.1) 

Gy(-) = y'(). 

System (4.1) can be written as: 

| Fr(-) = fCC) yC) 
Gy(-) =g, 2C), y()). 


Define an operator A: W?[0,7] >- WH[0, T], such that systems (3.1) and (3.2) 
can be represented as: 


(3.2) 


=> 


(3.3) AW =w, 

where A = [F,G]" and w = [zy]. 

Lemma 3.1. The operators F,G : W2[0,T] > W3[0,t] are bounded. 
Proof. Let x € W3[0,T]. Then 


2 
[Eel [0,T] -(Fe, Fa) aren 


(2.2) ator 
T 
=)? f (x""(t)) at 
E 2 
< (2(0))? + (2"(0))?+ f («""(t))°at 
0 


2 
= lal wW2[o,T]' 
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Similar argument can be used to show that the operator G from W2[0, T] into 
W3[0, 7] is bounded. 

Definition 3.2. Let j = 1,2. The space W’[0,7'] equipped with the inner product 

(RF 


? 


ko sae (£1, Y1)w [0,T] + (x2,2 dw: [0,7]? 


isa RK HS with reproducing kernel function Qas) = [a(-),¢s(-)]' where q,(-) 
is defined in (2.1) and R48) (4°) = [ri(-), rs()]" where r,,(-) is defined in (2.2), 
respectively. 


Lemma 3.3. The operator A: W7[0,T] > W'[0,T7] is bounded linear operator. 


Proof. Let X= [21,22]" belongs to W?[0,T] where x1, 22 € W3[0,T]. Then 


= 12 
JAX aon =|x ni 


=(x’, ies 


(24,04) iya tory t (22 bwapo.ry 
x 2 ri 2 g " 2 
=(2',(0))? + I (a(t) dt + (24(0))? + T (24) dt 
T 2 E 2. 
<(2r1(0))? + («4(0))? + i (2! (t)} dt + («2(0))? + (#4(0))? + f (wif (t))°at 
=lz1 lwo r]+lz2lw2to,r] 


AX enr 














Let {(£,7;)}#j-1 Þe a countable dense subset of D = [0,T] x [0,T] and define 
Vi = [hip] = FOE where A* is the adjoint operator of A and Bes is 
the kernel function of W'[0, T]. 


Lemma 3.4. Let {(&,7;)}?°_, be a countable dense set of points in D. Then Vi 


ijl 
is a complete system in W?[0,T], and 


Vij =ARG,,s;) 
for all i,j €N. 


Proof. Let Ù = [a,y]' € W?[0, 7] and let (T a a 0, for all i,j € N. 
Then ' 


(aT Wals, 1D) ya 0, T 


ha ,Qi;(s, a 
=AU (s;,tj) = 0, 


w![0,T] 
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where @(s;,tj;) = [x(si),y(t;)]". Since AW@(s;,t;) is vanish in a countable dense 
set of points in D, and it is continuous in W+[0, T] then w = 0. Furthermore, 


Vij =A" Q leir) 
-(A* Q (é:,7;) C, a) Ret,s;) C >) 


Nees JARg; aG Jha 
=AR¢,,5;)- 


W?[0,T] 


1[0,T] 





The proof is complete. 











In fact, Gram- Schmidt orthonormalization Process can be used to drive normal 
orthogonal basis {0; \2 of W?[0, T] from (Ware ©, as follows: 


Gi Pp 5: Bike, > Hat 


= T 


Bik Vex = [8:8], 


moes 


k=1 


where Big are orthogonalization coefficients of Yp. 


Theorem 3.5. Let {(£i,7;)}7j-1 be a countable dense subset of points of D and 


let X € W?[0, T] be a solution of AX = X is defined in (8.8). Then the solution 
X has the following representation: 


(3.4) 





SS BAX (te) ®; 
i=1 k=1 


=F Bie X'(ty) 
j=l k=1 

Since X (th) = [x (te), y' (th) 

plete. 


Ms 


re [S (tk Ek Yk), I(tr tk, e)] , so the proof is com- 














4. NUMERICAL RESULTS 


Example [3] A Lotka Volterra model is described by the following system of 
differential equations: 


dz = ax — bry 
(4.1) 
dy = -cy + dxy 


subject to the initial conditions z(0) = 1 = y(0), where a,b,c, and d are positive 
parameters 
e a: represents the natural growth rate of the pray in the absence of preda- 
tors; 
e b: represents the effect of predation of the pray; 
e c: represents the natural death rate of the predator in the absence of pray; 
e d: represents the efficiency and propagation rate of the predator in the 
present of pray. 
In numerical results we consider a = b = c= d= 1 and the implicit solution of the 
model is xy = exp(x + y - 2) 


5. CONCLUSION 
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